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SUMMARY 

Use  of  the  Monte  Carlo  SEMC  chemistry  data  package  has  authenticated  its 
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steady  state  experimental  data.  The  effects  of  finite  delay  on  avalanche  were 
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overestimates  of  the  avalanche  up  to  a  factor  of  about  2;  the  precise  error 
depends  on  the  time  interval  of  interest  and  on  the  ratio  E/p.  Having  proven 
the  accuracy  of  the  SEMC  chemistry  data  set  for  the  regime  where  experimental 
data  is  available,  we  believe  that  it  is  important  to  determine  the  effects  of 
avalanche  in  the  regime  E/p  >  (E/p)^j,.  Here  no  steady  state  exists  and 
experimental  data  is  not  available.  For  finite  time  intervals,  however, 
substantial  avalanche  can  be  expected.  Evaluation  of  avalanche  effects  in 
this  regime  for  the  pulse  widths  of  interest  is  an  important  next  step. 

The  data  from  the  Monte  Carlo  model  reported  by  Bloomberg  and  Pine  (198iJ) 
were  displayed  in  a  foi  m  useful  for  analysts  who  wish  to  check  the  importance 
of  avalanche  in  their  particle  push  algorithms.  Contour  plots  for  avalanche 
e-folding  values  of  f  =  1  and  f  =  5  are  presented  on  the  E-p  plane  for  a  given 
characteristic  time.  Alternatively,  an  expression  for  the  e-folding  value  for 
particular  E,  p,  and  characteristic  time,  tp,  is  given  by  Equation  (5). 
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CONVERSION  TABLE 
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1  000  000  X  E  -6 
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pound -force  (lbs  avoi  rdupoi.s) 

4.448  222 

newton  (N) 
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coulomb  Ailog  ram 

(CAg) 

shake 

1.000  000  X  E  -8 

second  (s) 

slug 

1.459  390  X  E  +1 

kilogram  (kg) 
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♦The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Dq  =  1  event/s. 
•♦The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 
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SECTION  1 

INTRODUCTION 
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It  has  become  clear  that  techniques  used  by  the  SGEMP  community  to 
predict  responses  in  tests  simulating  X-rays  incidents  on  missiles-in-flight 
are  inadequate.  A  major  reason  for  this  lack  of  fidelity  would  seem  to  be 
associated  with  the  absence  of  important  air  chemistry  in  algorithms 
describing  secondary  electrons  as  well  as  energetic  phtotelectrons.  A 
demonstrably  satisfactory  air  chemistry  data  set  has  recently  been  collected 
and  used  to  compute  secondary  electron  transport  in  a  Monte  Carlo  algorithm 
SEMC  (Bloomberg  and  Pine,  I98tl).  In  this  report  we  apply  results  from 
Bloomberg  and  Pine  (198*<)  to  quantify  the  effects  of  secondary  electron 
collisional  ionization  (avalanche)  over  the  short  time  intervals  of  interest 
in  the  study  of  X-ray  pulse  phenomena. 
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It  is  interesting  to  consider  why  the  community  has  so  long  eschewed  the 
incorporation  of  systematic  air  chemistry  data  into  the  algorithms  used  to 
predict  responses  for  missile-in-flight  lest  simulations.  The  reason  is  that 
complete  air  chemistry  data  is  relatively  unimportant  when  describing 
secondary  effects  in  the  very  low  and  very  high  pressure  regimes.  In  the  low 
pressure  regimes  secondaries  are  accelerated  rapidly  to  high  energies.  At 
high  energies  tne  electron  neutral  particle  interactions  are  relatively 
simple.  Large  angle  elastic  scattering  is  unimportant  and  the  energy  losses 
to  inelastic  processes  are  small.  On  the  other  hand,  at  high  pressure  the 
associated  air  chemistry  is  complicated.  Since  electron  distributions  quickly 
establish  steady  states,  however,  details  of  the  air  chemistry  are  not 
important  in  the  determination  of  electron  transport.  Adequate  empirical 
swarm  data  exist  in  this  regime  to  characterize  the  required  electron 
transport. 

It  is  clear  from  the  above  discussion  that  detailed  knowledge  of  air 
chemistry  is  required  in  order  to  characterize  the  intermediate  pressure 
regime.  The  implied  hope  of  the  community  that  the  effective  intermediate 
pressure  regime  would  be  sufficiently  narrow  so  as  not  to  affect  the 


phenomenology  seems  not  to  have  been  realized.  Indeed,  it  is  acknowledged 
that  important  phenomena  may  occur  just  in  the  intermediate  pressure  regime. 
Accordingly,  many  experiments  have  been  conducted  at  pressures  wherre  existing 
models  are  known  to  be  inaccurate.  The  work  of  Bloomberg  and  Pine  (198l^) 
incorporated  a  thorough  chemistry  data  set  for  nitrogen.  A  Monte  Carlo  model 
was  used  to  obtain  the  time  dependent  electron  distribution.  The  transient 
transport  coefficients  were  obtained  by  appropriate  particle  scoring.  The 
problem  was  restricted  to  the  case  of  spatially  uniform  distributions. 
Moreover,  the  electric  field  was  a  specified  fixed  value  throughout  space  and 
time.  Under  these  conditions  the  model  allows  a  consistent  evaluation  of 
electron  transport  coefficients  over  the  entire  range  of  pressures. 

This  report  begins  to  delineate  the  characteristics  of  the  important 
intermediate  pressure  regime  for  X-ray  pulses  incident  on  missiles-in- 
flight.  We  concentrate  our  attention  on  the  ionization  coefficient  or, 
alternatively,  the  collisional  ionization  frequency  for  secondary  electrons. 

In  Section  2  we  distinqulsh  the  cases  vSiere  secondaries  in  a  constant  electric 
field  can  either  approach  a  steady  state  or  a  runaway  condition.  In  Section  3 
we  examine  the  regime  for  which  steady-state  transport  is  possible.  We  argue 
that  there  is  a  delay  time  required  to  set  up  the  steady-state  distribution. 
For  finite  X-ray  pulse  durations  this  delay  can  change  the  avalanche  effect  as 
compared  to  when  steady  values  are  used  over  the  entire  pulse.  The  effects 
are  put  in  quantitive  terms.  In  Section  4  we  discuss  the  ramifications  of  the 
work  on  X-ray  simulations  of  missile-in-flight  environments.  In  particular, 
we  argue  that  the  effective  pulse  width  is  an  important  parameter  in  the 
representation  of  the  intermediate  pressure  regime.  The  results  are  discussed 
and  summarized  in  Sections  5  and  6,  respectively.  Appendix  A  contains 
mathematical  details  used  to  obtain  the  results;  Appendix  B  contains  the 
background  paper  of  Bloomberg  and  Pine  (1984). 
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SECTION  2 

STEADY-STATE  TRANSPORT  AND  SECONDARY  RUNAWAY 

It  is  important  for  us  to  characterize  the  division  between  the  high 
pressure  regime,  where  steady  state  obtains  and  the  low  pressure  regime,  where 
the  electric  field  continuously  accelerates  the  secondary  electrons.  The 
equation  of  motion  for  a  single  electron  is  given  by 

•"e  ^  ^  ‘"e^’m^  =  (1) 

Where  Vjj,  is  the  momentum  transfer  frequency,  that  is  the  frequency  for  large 
angle  elastic  scatter 


“  o„vnb  (2) 

where  is  the  energy  dependent  momentum  transfer  cross  section,  v  is  the 
electron  velocity,  and  rib  is  the  density  of  the  background  neutral 
particles.  The  largest  possible  value  of  the  drag  term  in  Equation  (1)  is 

">e"b  Maxfo^v"} 

It  is  reasonable  to  assume  that  whenever  the  electric  field  exceeds  this 
maximum  drag  that  electrons  will  run  away,  i.e.,  continually  gain  more  energy 
from  the  electric  field  than  is  lost  through  collisions.  The  critical  field 
so  defined  is  directly  proportional  to  the  background  number  density.  It  is 
convenient,  therefore,  to  define  a  critical  ratio  of  the  field  strength  to  the 
background  number  density  or,  alternatively,  the  field  to  pressure  ratio.  For 
an  N2  background  the  computed  critical  ratio  is 

(E/p)qp  =  1  .i|*10*V/(m-Torr)  (3) 
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Sworm  data  traditionally  are  plotted  against  E/p  as  the  independent 
parameter.  Since  the  data  are  steady  state  it  follows  that  swarm  theory  holds 
only  for  those  parameters  smaller  than  the  critical  value  of  Equation  (3). 
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SECTION  3 

TRANSIENT  BEHAVIOR 

The  feuiding  principle  in  the  work  of  Bloomberg  and  Pine  (198^1)  is  to 
include  sufficient  chemistry  de-excitation  as  well  as  elastic  and  ccllisional 
ionization  processes,  so  that  experimental  swarm  data  could  be  approximated  in 
the  time  asymptotic  limit.  A  measure  of  the  success  of  this  program  is  shown 
in  Figure  1 .  Here  the  ionization  rates  were  scored  over  the  time  histories 
tracked  in  the  Monte  Carlo  proceduie.  The  X*s  show  the  numerical  results  in 
the  long  time  limits  for  various  values  of  the  parameter  E/p.  The  solid  line 
in  Figure  1  shows  ionization  coefficients  from  the  experimental  swarm  data 
collected  by  Dutton  (1975).  The  agreement  with  experiment  is  good,  and  this 
is  verification  that  the  chemistry  incorporated  in  the  Monte  Carlo  algorithm 
is  sufficiently  complete  and  accurate  to  treat  the  time-dependent  electron 
transport  problem.  There  is  a  discrepancy  between  the  calculated  and  experi- 

S  _1  _1 

mentally  determined  ionization  cc  . 'icients  at  E/p  ~  7*10  V‘m  •  Torr  . 
This  may  indicate  that  the  critical  ratio  of  E/p  is  somewhat  lower  than  we 
have  estimated,  and  that  steady  state  is  not  being  achieved  at  the  upper  end 
of  the  curve  in  Figure  1 .  It  must  be  emphasized  that  the  agreement  between 
calculation  and  experiment  in  Figure  1  gives  us  confidence  in  using  the  Monte 
Carlo  algorithm  over  all  values  of  E/p.  For  values  above  the  critical  ratio 
there  will  be  no  steady  state,  but  the  resulting  transient  solutions  should 
give  the  correct  behavior. 

We  now  direct  our  attention  to  the  transient  behavior  of  the  frequency  of 
collisional  ionization  due  to  secondaries  for  values  of  E/p  <  (E/p)^,^,  i.e-., 
when  steady  state  is  achieved.  At  their  moment  of  creation  secondary 
electrons  (the  so-called  source  electrons)  have  low  energies.  In  the  presence 
of  an  electric  field  source  electrons  are  first  accelerated.  Eventually  the 
increase  in  electron  energy  from  the  field  is  balanced  by  losses  through 
collisions  with  the  background  neutral  particles.  Similarly,  a  momentum 
balance  is  achieved.  Obviously  there  is  a  delay  in  achieving  the  steady  state 
values.  This  is  shown  in  Figure  2,  where  the  ionization  coefficient,  a,  a 


measure  of  the  spatial  avaD^anohe  growth,  and  the  drift  velocity  are  plotted  as 
functions  of  time.  Note  that  the  delay  for  the  avalanche  coefficient  is  much 
greater  than  for  the  drift  velocity.  This  is  to  be  expected  since  the  drift 
velocity  is  a  measure  of  the  bulk  electron  distribution,  whereas  the 
ionization  coefficient  depends  on  the  distribution  of  electrons  with  energies 
large  enough  to  cause  ionization  {>15.6  eV  in  N2).  It  clearly  takes  longer 
for  the  wings  of  the  distribution  to  establish  their  steady  state. 

Collisional  ionization  by  secondaries  is  one  of  the  most  important 
processes  for  the  understanding  of  missile-in-flight  phenomena.  The  finite 
time  delay  in  reaching  steady  state  means  that  use  of  the  steady  stae 
ionization  frequency  over  the  entire  pulse  width  overestimates  the 
avalanche.  To  determine  the  effect  of  delay  we  prescribe  the  time  dependence 
of  the  ionization  frequency  by  function 

V  =  V  [1  -  exp  (-t/t^)3  (4) 

where  t^  is  the  decay  time  determined  from  best  fits  to  the  Monte  Carlo 
results;  v  is  the  steady  state  value.  We  show  in  Appendix  A  that  the  product 
pt^  is  itself  a  function  only  of  E/p.  The  functional  dependence  has  been 
determine  .  from  Monte  Carlo  runs  where  the  value  of  t^j  was  obtained  from  fits 
to  the  numerical  data.  The  resulting  plot  is  shown  in  Figure  3*  Note  that 
for  a  given  p,  t^  decreases  with  increasing  E.  This  result  is  expected  since 
the  higher  field  acceleration  should  diminish  the  time  required  for  initially 
stationary  secondary  electrons  to  achieve  their  steady-state  distribution. 

The  pressure  normalized  steady-state  collision  frequency  is  also  only  a 
function  of  E/p.  The  functional  form  can  be  recovered  from  the  Dutton  (1975) 
data  and  is  shown  in  Figure  4.  This  completes  specification  of  the  time 
dependent  ionization  frequency  for  values  of  E/p  small  enough  so  that  a  steady 
state  value  eventually  is  achieved,  as  implied  in  Equation  (4). 


SECTION 


RESULTS 

We  now  discuss  the  effects  of  non-zero  time  delays  over  finite  pulse 
widths  tjj.  All  thin  air  SGEMP  processes  are  limited  to  finite  time 
intervals.  The  interval  of  concern  here  is  that  for  which  the  electric  field 
strength  remains  relatively  constant.  This  duration  is  difficult  to 
determine,  because  it  is  affected  in  part  by  self-consistent  effects  from  the 
very  ionization  phenomena  we  are  attempting  to  describe.  Nevertheless,  it  is 
useful  in  this  analysis  to  introduce  an  effective  time  interval.  In  general, 
the  time  of  interest  can  be  associated  simply  with  the  X-ray  radiation  pulse 
width. 

The  importance  of  incorporating  finite  delay  time  into  calculations 
depends  on  the  time  interval  t^j.  It  is  obvious  from  Equation  (li)  that  neglect 
of  the  delay  time  overestimates  c.he  number  of  ionizations  in  any  given  time 
interval.  The  magnitude  of  the  overestimating  error  resutling  from  use  only 
of  the  steady-state  ionization  frequency  over  the  time  interval  has  been 
calculated.  Details  of  the  calculation  are  presented  in  Appendix  A,  and  the 
results  are  plotted  in  Figure  5.  Here  the  overestimate  factor  for  several 
values  of  the  parameter  ptp  is  plotted  as  a  function  of  E/p.  For  small  ptp 
the  overestimate  factor  is  close  to  unity  over  the  range  of  E/p  under 
consideration.  This  is  simply  a  manifestation  of  the  fact  that  the  time 
interval  is  so  small  that  hardly  any  collisional  ionizations  result  from  a 
secondary  electron.  Thus  the  delay  time  has  little  effect,  since  ionization 
is  unimportant  over  the  time  interval  even  when  the  steady-state  value  is 
used.  For  larger  ptp  neglect  of  the  delay  time  leads  to  a  more  significant 
error.  Finally,  for  very  large  ptp  the  overestimate  factor  reaches  a  limit 
for  a  given  value  of  E/p.  This  limit  is  easily  understood.  For  t>>t^  the 
actual  ionization  rate  is  nearly  the  steady-state  value.  Thus,  increasing  the 
value  of  tp  when  tp>>t^  will  not  lead  to  additional  increases  in  the 
overestimate. 
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It  is  interesting  to  note  that  the  overestimate  does  not  exceed  a  value 
of  about  2,  which  occurs  for  large  E/p.  The  reason  for  this  is  that  while  the 
ionization  frequency  becomes  important  for  large  E/p,  the  delay  time  decreases 
for  large  E/p,  so  that  the  importance  of  delay  is  reduced.  The  net  effect  is 
that  the  overestimate  is  limited  to  a  rather  modest  value.  A  more  detailed 
discussion  of  this  result  is  given  in  Appendix  A. 

Having  shown  that  consideration  of  the  delay  time  is  important  in  the 
calculation  of  the  number  of  secondary  collisional  ionization  events  over  a 
pulse  interval.  Ip,  we  now  address  the  larger  question  of  when  secondary 
collisonal  ionization  is  an  important  effect.  Collisional  ionization  is  an 
important  effect  whenever  the  corresponding  e-fold  vdt  =  r  ^  1 ,  where  v  is 
given  by  Equation  (J4). 

In  Appendix  A  we  show  that  for  any  selected  e~fold  r  a  curve  in  E/p  -  ptp 
space  is  defined.  In  Figure  6  we  present  two  curves  for  both  a  modest  e-fold 
(r=1 )  and  a  large  e-fold  (r=5)  for  collisional  ionization.  In  order  to  get  a 
particular  avalanche  multiplication  over  the  time  interval  for  smaller , ptp  the 
value  of  E/p  must  increase.  For  convenience  consider  p  fixed;  then  as  t 
decreases  a  larger  avalanche  rate  is  required  to  give  the  same 
multiplication.  As  E  gets  to  be  too  large  a  runav/ay  condition  applies,  and  we 
no  longer  can  apply  the  present  theory.  On  the  other  hand,  for  very  large  tp 
a  smaller  avalanche  rate  is  required  for  the  same  multiplication.  The  value 
of  E  does  not  decrease  indefinitely,  however,  since  there  is  a  critical  value 
of  E/p  below  which  there  is  a  zero  avalanche  rate. 

We  have  unfolded  the  results  of  Figure  6  in  Figure  Ta-j  to  give  e-fold 
contour  plots  for  r=1  and  r=5  onto  the  E-plane  for  fixed  values  of  tp.  This 
format  is  a  useful  one  to  check  the  validity  of  the  treatment  of  secondary 
collisional  ionization  in  the  numerical  algorithms  used  for  this  air  SGEMP 
predictions.  For  these- predictions  the  pressure  is  specified.  The  codes 
themselves  indicate  the  E  field  strengths  and  durations.  Thus,  the  adequacy 
of  air  chemistry  data  in  the  predictive  codes  can  be  checked  straight¬ 
forwardly. 
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Let  us  now  discuss  the  contour  plots  in  the  E-plane  when  we  fix  the  pulse 
tjj.  Note  as  we  increase  p  the  required  value  of  E  increases.  The  universal 
plot  of  Figure  6  indicates  that  on  a  given  e-fold  contour  the  ratio  E/p 
decreases  for  increasing  p.  This  implies  that  the  increase  in  the  required  E 
with  p  is  sub-linear,  and  this  conforms  with  the  plots  on  Figures  7a-j.  For 
large  enough  p  the  required  value  of  E  increases  most  rapidly  (nearly  linearly 
with  p) .  This  occurs  because,  while  the  avalanche  rate  required  for  a  given 
multiplication  does  not  change,  the  I'atio  v/p  decreases  with  increasing  p. 

For  small  enough  v/p  the  corresponding  value  E/p  is  essentially  constant,  so 
that  E  increases  nearly  linearly  with  increasing  p. 

In  summary.  Figures  7a-j  give  an  indication  (for  a  given  characteristic 
time  and  pressure)  the  extent  to  which  a  given  E-field  can  lead  to  significant 
collisional  ionization  by  secondary  electrons.  It  may  be  useful  to  compare 
quantitatively  the  secondary  collisional  ionization  from  self-consistent  thin 
air  SGEMP  predictive  codes  with  the  results  of  this  work.  Fc*’  this  reason  we 
reference  (from  Appendix  A)  the  expressions  for  the  growth  facu^.  of  a  single 
initiatory  secondary  electron. 

The  growth  or  multiplication  factor  is  G  =  exp(r)  where  r  is  the  e-fold 
defined  by 

P  =  g(E/p)  Cptp  +  f(E/p)exp{-ptp/f (E/p)}  -  f(E/p)]  (5a) 

where 

f(E/p)  =  pt^  (5b) 

and 

g(E/p)  =  v/p  (5c) 

are  plotted  in  Figures  3  and  respectively. 


SECTION  5 


DISCUSSION  OF  RESULTS 

Determination  of  the  importance  of  avalanche  in  thin-air  SGEMP  requires 
information  about  the  electric  field  strength  throughout  the  air  around 
missile  surfaces.  The  electric  field  has  not  been  a  parameter  of  physical 
concern  in  most  past  analyses  of  thin-air  SGEMP  phenomena.  Rather  a  greater 
emphasis  has  been  placed  on  the  determination  through  theory  and  experiment  of 
the  voltages  at  various  missile  surfaces.  As  mentioned  before,  the  emphasis 
on  differential  voltages  has  been  motivated  by  the  belief  in  the  community 
that  the  net  flow  of  primary  photoelectrons  from  one  section  of  the  missile  to 
another  is  the  most  important  thin-air  phenomenon.  In  addition  to  stressing 
the  importance  of  the  electric  field,  we  have  also  stressed  the  role  of  the 
characteristic  time  interval  for  the  field.  It  is  this  quantity  (or  more 
generally  the  product,  pt^)  that  determines  the  importance  of  the  delay  time 
in  the  computation  of  the  collisonal  ionization  frequency. 

The  issue  of  the  duration  of  the  effective  time  interval  for  the  electric 
field  is  complex  but  important.  There  are  two  aspects  to  the  problem.  The 
first  is  structural  and  concerns  how  various  parts  of  a  modeled  system  connect 
with  each  other  electrically.  The  second  concerns  air  itself  where  there  is 
another  relaxation  channel  due  to  the  number  of  secondary  charge  carriers. 

This  is  determined  by  the  number  of  initiating  electrons  as  well  as  by 
avalanche  itself.  In  air  directly  exposed  to  incident  X-rays  the  creation  of 
electrons  through  photoemission  or  by  collisional  ionization  of  photoelectrons 
emitted  from  surfaces  may  be  very  large.  In  such  cases  the  air  conductivity 
might  be  sufficiently  large  so  as  to  rapidly  reduce  field  strengths  through 
relaxation.  Then  there  might  be  insufficient  time  for  significant  avalanche 
to  occur.  Another  possibility  is  that  only  a  moderate  number  of  initiating 
secondary  electrons  are  produced  in  air.  This  could  happen  if  shadowing  is 
important.  The  conductiv.',ty  due  to  the  initiating  electrons  would  be 
relatively  small  so  that  the  field  strengths  could  be  maintained  for  long 
times.  Under  these  circumstances  avalanche  could  become  important. 


The  present  work  extends  naturally  to  include  causes  for  E/p  >  (E/p)^^. 

In  this  regime  there  is  no  steady  state,  and,  after  long  times,  the  collisonal 
ionization  frequency  decreases  monotonically  with  increasing  time.  Over  the 
finite  time  intervals  of  interest,  however,  the  total  avalanche  for  E/p  > 
(E/p)qj,  may  be  substantial. 
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solid  curve  shows  experimental  data  assembled 
corresponding  values  obtained  in  this  work  ar 


Figure  2.  Drift  velocity  and  coefficient  as  functions  of  time 


Lotted  as  a  func 


steady  State  Overestimate  Factor 


These  plots  show  the  error  in  the  nmber  of  ionizations 
occurring  in  a  pulse  interval,  t  ,  in  seconds  when  the 
ionization  time  delay  is  neglectld. 
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Figure  5.  Steady-state  overestimate  factor. 
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APPENDIX  A 


The  scaling  relations  associated  with  the  transient  expression  for  the 
collisonal  ionization  frequency  caused  by  secondary  electrons  can  be  derived 
through  general  arguments.  The  equation  for  the  electron  distribution 
function  for  a  uniform  system  with  a  constant  electric  field  E  can  be  written 
symbolically  as 


§1  =  E‘G  -  N‘H  (A-1 ) 

where  N  is  the  neutral  number  density  and  G  and  H  depend  on  the  electron 
velocities  and  not  explicitly  on  time.  While  Equation  (A-1)  generally 
represents  a  Boltzmann  equation,  it  can  also  represent  the  simple  electron 
equation  of  motion  often  used  in  heuristic  descriptions  of  many  particle 
systems.  In  this  case 


F  ^  V 
G  •*  e/m 
H  v'v 


where  v  is  the  particle  velocity,  e/m  the  charge  to  mass  ratio  for  an 
electron,  and  v'  the  secondary  electron  collision  frequency  per  single  neutral 
particle  in  a  unit  volume;  the  conventional  drag  frequency  Vjjj  =  v*N. 

It  is  simple  to  put  Equation  (A-1)  in  the  form 

^  -  H  (A-2) 

where  C  =  Nt.  It  is  immediately  seen  from  Equation  (A-2)  that  for  a  given  E/N 
(or  equivalently  E/p,  where  p  is  pressure)  the  time  evolution  for  F  varies 
inversely  with  N  (or  p).  In  steady  state  F  becomes  independent  of  C  and 
depends  on  E/p  alone. 


The  transient  collisional  ionization  frequency  has  been  calculated  by 
Bloomberg  and  Pine  (1984).  Time  variation  of  the  collision  frequency  can  be 
approximated  by  the  expression 

V  =  V  (1  -  (A-3) 

where  t^  is  a  delay  time  obtained  from  the  numerical  results  of  the  Monte 
Carlo  algorithm;  \5  is  the  steady  state  collisional  ionization  frequency.  From 
the  scaling  results  discussed  above  we  can  write 

t^  =  f{E/p)/p  (A-4) 

V  =  p  g(E/p)  (A-5) 

The  collisonal  ionization  e-folding  over  a  time  interval,  tp,  is 

r  =  /  P  vdt  (A-6) 

o 

Substitution  of  Equation  (A-3)  into  Equation  (A-6)  and  use  of  Equations  (A-4) 
and  (A-5)  gives 


r  =  8  [pfcp  +  f  exp{-ptp/f}  -  f]  (A-7) 

which  expression  was  reproduced  in  the  main  Text  as  Equation  (5a). 

The  increase  in  the  number  of  electrons  created  in  the  avalanche  of  a  single 
secondary  electron  over  a  time  interval  tp  is  just  exp(r).  Use  of  the  steady 
state  value  for  the  collisional  ionization  over  the  entire  time  interval  tp 
overestimates  the  avalanche  by  the  factor 

exp{gptn  -  1’}  =  exp{fg  [1  -  exp(-pt„/f )]}  (A-8) 
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Abstract 

A  Honte  Carlo  algorithm  Is  used  to  determine 
the  time  behavior  of  source  secondary  electrons  for 
ranges  of  the  electric  field  to  pressure  ratio  £/p  of 
interest  In  AHSGEHP.  The  algorithm  contains  a  very 
detailed  cross  section  set  describing  electron 
interactions  with  the  bacltgraund  gas.  We  show  that 
the  delay  in  the  attainment  of  the  peak  tine 
independent  ioniaation  frequency  (or  ionization 
coefficient)  nay  result  in  negligible  loniaaticn  over 
times  of  interest.  In  any  case  the  behavior  is  shown 
to  behave  much  differently  than  in  examples  where 
limited  cross  section  sets,  common  in  currently 
employed  predictive  codes,  are  employed.  In 
particular,  the  importance  of  momentum  transfer  is 
indicated.  A  critique  of  the  scaling  inplications  of 
the  phenomena  is  made. 


1.  Introduction 


The  study  of  air  modified  SGEHP-.AKSGEHP 
requires  the  understanding  of  secondary  electron 
behavior  over  a  wide  range  of  electric  fields  and 

ambient  pressure^'^.  The  maximum  electric  field, 
associated  with  the  Chi1d>Langmuir  space  charge  limit 

for  photoemitted  electrons.  Is  about  5  x  10^  V/m, 
whereas  the  pressure  range  for  which  X>ray  penetration 
is  important  is  vacuum  to  about  10  Torr.  These  numbers 
imply  a  wide  range  for  the  electric  field  to  pressure 
ratio  E/p,  the  traditional  parameter  against  which 
electron  swarm  transport  quantities  are  specified.  A 
last  general  parameter  of  interest  in  AHSGEHP  problem 
is  the  radiation  pulse  width,  which  ranges  through 
tens  of  nanoseconds. 

While  the  behavior  of  secondary  electrons  is 
essential  to  AHSGEHP,  the  models  developed  thus  far 
are  not  adequate  to  describe  these  electrons  for  the 
wide  range  of  parameters  of  interest.  Indeed,  the 
models  describe  accurately  only  two  extremes  in 
parameter  space.  The  first  is  for  low  E/p,  where  the 
secondaries  approach  a  steady  state  distribution 
withi.'i  times  much  shorter  than  the  radiation  pulse 

width.  In  this  case  the  swarm  theory^  description  is 
applicable.  The  second  extreme  is  for  high  E/p  where 
each  secondary  reaches  a  very  large  energy  and  all 
interactions  except  for  collisional  ionization  can  be 
neglected.  In  this  case  no  distinction  is  made 
between  secondary  and  high  energy  primary  photoemitted 
electrons. 

In  this  work  we  determine  the  transient 
properties  of  the  secondary  electron  population  caused 
by  initiatory  electrons  moving  within  given  constant 
electric  fields  and  interacting  with  N2  at  arbitrary 

pressures.  The  problem  is  addressed  through  a  Monte 
Carlo  approach.  Honte  Carlo  techniques  have  been  used 
extensively  to  study  the  transient  effects  of  electron 
Swarms  in  the  pre-breakdown  phase  of  gas 
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discharges  .  Besides  applying  Honte  Carlo  to 
parameters  of  interest  in  AHSGEHP,  we  have  generalized 
previous  work  in  several  ways.  First,  we  have 


incorporated  into  the  algorithm  accurate  cross  section 
sets  for  energy  loss  due  to  excitation  as  well  as, 
collisional  ionization.  In  addition,  we  include  the 
effects  of  elastic  scattering  through  use  of  momentum 
transfer  cross  sections.  We  also  have  speeded  up  the 
previous  calculations  for  determination  of  the  mean 
free  tiroes  for  a  particle  history,  when  the  particle 
energy  changes  appreciably  between  collisons. 

Finally,  we  have  employed  sampling  techniques  to 
enable  us  to  find  important  history  averaged 
Information  such  as  velocity  drift  and  the  ionization 
coefficient  with  the  use  of  a  relatively  small  number, 
of  histories.  Thus,  useful  cvvnputations  were  carried 
cut  with  nominal  computer  times  even  with  the  detailed 
cross  section  sets. 

The  algorithm  is  described  in  Section  2.  In 
Section  3  it  is  shown  how  the  secondary  temperatures, 
drift  velocity  evolve  to  steady  state  as  a  function  of 
electric  field  and  pressure.  We  benchmark  the  steady 
state  values  for  the  ionization  coefficients  against 

c 

the  tabulated  results  of  Dutton  for  various- E/p.  We 
show  how  the  approach  to  steady  state  depends  on  the 
electric  field  even  though  steady  state  values  depend 
only  on  E/p.  In  Section  4  the  results  of  tlie  Honte 
Carlo  algorithm  are  shown  when  a  limited  cross  section 
set,  characteristic  of  the  interactions  employed  in 
current  A.HSGEKP  computer  simulations,  is  used  in  place 
of  the  complete  set.  Results  for  the  two  cross 
sections  are  compared.  In  Section  5  results  are 
discussed.  We  define  the  limits  of  applicability  of 
A.HSGEHP  algorithms  as  a  function  of  the  X-ray  pulse 
width.  We  indicate  the  importance  of  the  errors 
resulting  from  the  limited  cross  section  set. 


2.  Single  Scatter  Honte  Carlo 


The  algorithm  described  here  was 

incorporated  into  an  existing  code  SEHC^^  (Single 
Electron  Honte  Carlo)  developed  to  solve  for  electron 
properties  in  the  limit  of  large  enough  times  such 
that  transport  coefficients  can  be  defined.  Here  «e 
discuss  those  changes  needed  to  compute  time  dependent 
properties  of  secondary  electrons  in  Nj.  Hany 

features  of  SEHC  were  retained  in  the  modification  for 
explicitly  transient  behavior.  Among  these  were  the 
scorings  to  determine  the  velocity  distribution 
functions,  as  well  as  the  scattering  kernel 
formulation.  The  last  procedure  connects  the  Honte 
Carlo  approach  to  the  Boltzmann  equation.  In  a 
practical  sense,  it  also  is  given  a  convenient  format 
for  the  input  specification  of  detailed  tabulated 
cross  Section  sets  for  arbitrary  background  materials. 

Algorithm  for  Transient  Effects. 

The  proble.n  is  to  determine  secondary 
electron  prooerties  evolving  fran  a  set  of  initiatory 
electrons  existing  within  an  background  in  the 

presence  of  a  constant  electric  field.  In  the  spirit 
of  the  single  scatter  Honte  Carlo  approach  we  can 
imagine  a  single  electron  being  launched  at  time  zero. 
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The  electron  undergoes  cullUlons  wUh  the  background 
during  the  course  of  tts  history  ending  after  a 
specified  tiiae  interval.  For  a  sufficiently  large 
energy  the  electron  creates  an  ion  pair  through 
coilisional  ionization.  This  second  particle  created 
some  time  after  time  zero,  should  also  be  tracked  by 
the  algorithm.  Since  all  histories  end  at  the  same 
time,  the  lifetime  of  the  second  electron  is  shorter 
than  the  first.  More  particles  are  created  as  the 
avalanche  proceeds.  Within  the  time  interval  of 
interest  a  finite  number  of  electrons  will  be  created, 
and,  in  principle,  all  of  them  could  be  tracked.  To 
improve  statistics  a  second  electron  could  be  launched 
at  time  zero,  and  the  entire  process  could  be 
repeated.  The  difficulty  with  this  straightforward 
selection  of  electron  histories  is  that  the  total 
number  of  histories  '■annot  be  chosen  prior  to  the 
calculation.  The  cost  of  such  a  run  cannot  be 
controlled  a  priori.  To  control  this  problem  we 
specify  initially  the  number  of  time  histories  to  be 
tracked  in  a  code  run.  The  first  particle  history  is 
launched  at  time  zero.  We  count  the  number  of 
ionizations  created  by  the  particles  and  thereby 
obtain  a  first  estimate  to  the  ionization  rate,  time 
averaged  over  the  time  interval  of  interest.  The 
second  particle  history  commences  at  either  time  zero 
or  at  one  of  the  times  when  an  ionization  event  took 
place.  The  choice  for  the  launch  time  of  the  second 
particle  Is  determined  by  a  Monte  Carlo  selection  from 
the  distribution  of  particles  over  the  complete  time 
interval  determined  by  avalanche  statistics.  This 
means  that  electrons  created  by  coilisional  ionization 
late  in  time  have  a  higher  probability  of  being 
selected  for  a  launch.  The  procedure  is  then 
repeated.  The  initial  times  for  candidate  particle 
histories  are  time  zero  and  those  limes  for  which 
ionizations  occurred  in  all  previous  histories.  The 
distribution  used  to  determine  the  time  for  each 
launch  is  determined  from  statisticil  information 
obtained  in  all  previous  histories.  For  simplicity, 
this  analysis  considers  that,  at  the  beginning  of  each 
time  history,  the  electron  is  at  rest. 

For  sufficiently  high  E/p  a  number  of 
particles  may  run  away  over  the  time  interval  of 
interest.  Runaways  are  a  component  of  the  secondary 
electron  population  that  are  difficult  to  incorporate 
into  simple  swarm  theory.  In  large  constant  electric 
fields  they  continually  take  energy  from  the  field. 

In  a  realistic  situation  the  energy  that  the  electric 
field  can  give  to  an  electron  must  be  limited,  but  the 
assignment  of  a  limit  necessarily  brings  into  play  the 
parameter  of  a  particular  configuration.  In  AHSGEMP 
problems  we  have  selected  10  keV  as  the  limiting 
secondary  electron  energy.  This  corresponds  to  about 
the  largest  potential  difference  that  can  be  sustained 
by  the  primary  photocmitted  electrons.  For  our 
purposes  the  transient  swarm  properties  are  not 
sensitive  to  the  prescribed  runaway  cut  off  energy, 
since  multi  keV  electrons  contribute  negligibly  to  the 
ionization  coefficients  when  the  bulk  of  the 
secondaries  have  energies  less  than  1  keV. 

In  Figure  1  we  show  a  scnematic  diagram 
Illustrating  the  types  of  time  histories  considered 
over  a  specified  time  interval  t^.  Only  those 

particles  that  .end  their  histories  at  t^  contribute  to 

the  velocity  distribution  function,  that  is  no 
velocity  information  from  runaways  is  retained.  On 
the  other  hand,  the  contribution  to  coilisional 
ionization  of  all  histories  prior  to  possible  runaway 
are  included. 

Collision  Cross  Sections. 

Elastic  and  inelastic  collision  cross  sections 


have  been  incorporated  into  the  numerical  algorithm  to 
describe  electron  interaction  with  the  background 

Besides  ionization  important  excitation  cross  sections 
were  included  in  the  set  of  inelastic  cross  sections. 
These  data  were  obtained  from  sources  compiled  by 

Strickland  et  al.^^ 
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Figure  1.  Four  Types  of  Electron  Histories 


A.  An  initiatory  electron  launched  at  t  c  0  and 
ending  its  history  at  the  prescribed  time  t  »  t^; 

8.  An  electron  created  at  t  0  in  the  avalanche  and 
ending  its  history  at  t  «  t^; 

C.  An  initiatory  electron  that  runs  away  at  t  <  t^; 

D.  An  avalanche  created  electron  running  away  at 


The  momentum  transfer  coefficients  were 
computed  from  differential  elastic  cross  section  sets 

from  Siyn  et  al.  below  70  eV,  and  from  Bromberg 
below  500  eV.  The  momentum  transfer  cross  section  is 
rel.'.ted  to  the  differential  cross  section  by  tne 
formula 


Om(0)  =  2ii  ■ 

ni'  ' 


In  general  a  numerical  computation  of  this  equation 
must  be  carried  out.  For  energies  above  1  keV  we  use 
the  form  for  the  screened  Rutherford  potential  used  by 

Strickland  et  al.^^. 

The  cross  sections  were  used  to  determine 
coilisional  events  for  each  electron  history.  The 
particle  collision  frequenc"  v(U)  is  related  to  the 
sura  of  all  cross  sections  'J)  for  a  particular  energy 

U  by  the  expression 


/S 


^  Siny(2  -  cosv)  (J*. 


(1) 


v.f(U)  ■=  Oy(U)  •  Vau/m  (2 

The  time  for  a  collision  to  occur  is  related  to  a 
random  number  c  between  0  and  1  and  me  mean  frequency 
by  the  formula 

ln(l  -  0  =  -  J 

0 
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where  v^  ts  evtlutted  for  the  chtnging  energy  of  the 

electron  experiencing  between  collisions  the  constant 
acceleration  from  the  electric  field  directed  along 
the  r-exis,  E.  Once  the  logarithm  is  determined  by 
selection  of  the  random  number  the  integration  is 
carried  out  in  time  until  the  prescribed  value  of  the 
integral  Is  reached.  This  procedure  defines  the  time 
interval  until  the  next  collision.  It  is  relatively 
easy  to  ^rform  even  when  vt  varies  considerably  over 
the  integration  range.  ' 

Once  the  time  fnr  the  next  collision  has  been 
established  one  needs  to  determine  the  type  of 
collision.  This  is  obtained  by  consideration  of  the 
relative  weignt  of  all  cross  sections.  The  specific 
collision  type  is  found  by  random  number  selection. 

The  effect  of  inelastic  collisions  is  that  the 
particle  loses  a  specified  energy  with  no  directional 
change  of  velocity  vector  while  particles  undergo  a 
change  of  direction  based  upon  an  isotropic 
distribution  when  tney  transfer  Romenttn. 

In  Figure  2  we  sunnarize  the  cross  section 
information  as  a  function  of  electron  energy.  The 
pressure  nornalized  momentum  transfer  frequency  v„/p 

is  defined  in  terms  of  nooentum  transfer  cross  section 
o_  by  the  relation 


VP  ‘  ‘’m 


(4) 


where  ii  is  the  number  density  of  N,  corresponding  to  a 
-  IR  1 

pressure  of  1  torr:  N  »  3.55  x  10  cm'  ,  v  is  the 
electro.-;  velocity.  The  second  plot  is  the  energy  loss 
rate  normalized  to  pressure.  Here  v^/p  is  related  to 

the  loss  function  L(£)  by  the  relation 


Vy/P 


vl.(U)/U  •  N 


(5) 


where  the  loss  function  represents  the  energy  loss 
from  all  important  channels.  In  our  case  we  consider 
several  collisional  excitation  loss  mechanisms  as  well 
as  the  collisional  ionization  energy  loss  channel. 
Therefore, 


m 


(i) 


(6) 


where  the  sum  extends  over  all  important  excitation 
channels;  ij^s  the  cross  section  for  the  i^**  type  of 
excitation  while  is  the  energy  lost  through  this 
channel,  o,  is  the  ionization  cross  section  and  c,  is 


“I 

the  energy  lost  in  ionization,  thus  In  N21 
15.6  eV. 

3.  Solutions 


In  Figure  3  we  plot  the  cumulatively  averaged 
ionization  coefficient  and  the  instantaneous  drift 
velocity  as  functions  of  time.  We  have  selected  the 
steady  state  value  of  the  ionization  coefficient  as 
the  time  dependent  transport  coefficient  quoted  In  the 
electron  swarm  literatire.  We  have  compared  the 
steady  state 

loniz.tion  coefficients  for  various  value»  of  E/p  with 

5 

the  results  of  a  data  survey  conducted  by  Dutton  . 


The  good  agreement  shown  In  figure  4  gives  credence  to 
the  accuracy  of  the  detailed  cross  section  set  used  in 

the  analysis.  The  discrepancy  at  t/p  •  10^  V  x  m'*  x 

torr'^  is  not  uhexpecied.  First  it  occurs  at  a  high 
enough  value  of  E/p  where  runaways  might  occur  so  tnat 
the  data  will  be  less  reliable.  Secondly,  the  loss 
function  employed  in  the  analysis  neglects  the  kinetic 
energy  Imparted  to  the  electrons  created  by 
collisional  ionization.  This  source  of  energy  loss 
for  Ionizing  electrons  becomes  Important  for  high 
electron  energies. 


lectron  Energy 


Figure  3.  Drift  Velocity  and  Cumulatively  Averaged 

Ionization  Coefficient  as  Functions  of  Time. 


Parameters  are  E. 


2  X  lO^V  X  m'^ 


torr.  500  histories  were  launched. 


Although  the  asymptotic  transport 
coefficients  for  electron  swarms  depend  only  on  the 
parameter  E/p,  the  rate  at  which  these  coefficients 
reach  the  asymptotic  value  depends  on  the  value  of  the 
field.  In  fact,  consideration  of  the  Boltzmann 
equation  equivalent  to  the  Monte  Carlo  procedure  shows 
that  there  exists  a  simple  scaling  for  the  rise  times. 
Assume  that  the  rise  time  for  some  property  computed 
for  a  given  E  and  p  is  t  Then  for  the  case  where 
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another  field  v*1ue  E'  •  EE  *nd  another  pressure  p'  > 
kp,  Mhere  k  is  any  constant,  the  corresponding  rise 
tise'  is  related  to  the  computed  time  by  the  relation 
tp  ■  t|./k.  The  Honte  Carlo  algorithm  produces  the 

required  scaling.  This  is  shown  in  Figure  S  where  the 
pressure  normalized  ionization  frequency  is  plotted 

against  time  for  E  •  2  x  10^  V  x  and  E  °  2  x  10^ 

V  X  when  p  «  10'^  torr  and  1  torr,  respectively. 


Figure  4.  Steady  State  Ionization  Coefficient  vs  E/p 
(from  Dutton).  The  values  obtained  in  this 
work  are  shown  by  the  X's. 


rigure  5.  Cumulative  Ionization  Rate  vs.  Time  for  E  • 

2  X  )0^  and  E  •  2  X  iO'  V  X  m"^.  In  each 
4  -1 

case  E/p  «  2  X  10  X  m  X  tor>-  1, 


4.  Comparison  with  limited  Cross  Section  Set 


The  particle  pushing  algorithms  used  for 
experimental  predictions  in  the  AMSGEHP  cowiunity 
generally  exclude  the  effects  of  elastic  scattering 
and  energy  loss  through  discrete  excitation. 

Generally,  the  only  electron  background  interaction 
included  is  that  for  collisional  ionization.  To  test 
the  effect  of  excluding  the  elastic  and  excitation 
interactions  we  have  limited  the  original  cress 
section  set  to  include  only  the  collisional  ionization 
cross  section.  The  results  of  those  runs  were 
compared  to  runs  with  the  same  parameters  when  the 
complete  cross  section  set  was  employed. 

The  results  are  plotted  in  Figures  6  and  7. 

In  Figure  6  we  show  the  drift  velocities  as  functions 
of  time.  >'t  is  apparent  that  the  limited  cross 
section  set.  gives  rise  to  drift  velocities  that  reach 
steady  state  values  about  6  times  higher  than  the 
corresponding  drift  velocities  obtained  with  the 
cofflclcte  cross  section.  The  values  of  the  drift 
velocity  for  tne  latter  case  agrees  well  with  the  data 

reported  by  Dutton^.  A  large  discrepancy  is  also  shown 
in  Figure  6  where  the  cumulatively  averaged  ionization 
coefficient  is  plotted  against  time.  The  ionization 
coefficient  for  the  limitedjSet  peaks  early  at  a  much 
higher  value  than  for  the  case  of  the  complete  set. 

The  steady  state  values  achieved  for  this  case  again 
agrees  well  with  the  data  of  Dutton  (see  Figure  4). 


tlXr  (l) 


Figure  6.  Drift  Velocity  Versus  Time 

The  solid  curve  is  the  numerical  resuU  for  the 
complete  uderacticn  cress  section  set,  while  the 
dashed  curve  is  tne  result  for  the  limited  cress 

C  1 

section  set.  Parameters  are  E  »■  10  V  x  m  ,  p  »  10 


5.  Discussion  of  Results 


We  have  shown  large  discrepancies  in 
Ionization  coefficients  (and  lontzalion  rates)  between 
the  results  with  a  complete  cross  section  and  with  a 
limited  sat  approximating  the  treatment  of  secondaries 
in  present  AHSGEK?  predictions  codes.  Tiie 
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discrepinctes  tre  aost  tevere  it  1m  E/p  t^ere  the 
aonentbo  trinsfer  rites  ind  excUitton  processes  ire 
very  {aportint.  At  higher  E^,  where  the  secondirles 
hive  lirger  energies,  excltitlon  processes  ind 
Momentum  trinsfer  rites  no  longer  strongly  donlnite 
the  collision  lonlzitlon  processes.  Thus,  the 
dlscrepincies  in  the  results,  while  significint,  are 
not  IS  severe. 


Figure  7.  CurouUtively  Avenged  lonizition  Coefficent 
Versus  Tine. 

The  solid  curve  is  the  numerical  result  for  the 
complete  interaction  cross  section,  while  the  dashed 
curve  is  the  result  for  the  limited  cross  section 

set.  Parameters  are  E  «  10®  V  x  m'\  p  «  10  torr. 


The  fundamental  question  to  be  asked  in  .WS6EHP 
phenomena  is  whether  the  discrepancies  affect  the 
physics.  The  key  parameter  Is  the  ratio  of  the 
ionization  rate  with  the  characteristic  frequency  of 
the  radiation  pulse.  If  this  ratio  is  small,  than  the 
avalanche  may  not  be  an  important  effect,  since  the 
initiatory  electrons  are  being  generated  at  a  rate  of 
order  the  inverse  of  the  characteristic  radiation 
pulse  widiii.  Since  the  avalanche  ionization  rate  is 
small  for  1m  E/p,  the  large  discrepancies  observed 
when  one  uses  a  complete  cross  set  compared  to  the 
limited  set  nay  not  affect  the  physics.  This  is 
particularly  true  for  a  given  E/p  when  p  is  small. 

Then  the  transient  rise  time  of  the  ionization  rate  to 
its  asympotic  value  may  be  even  larger  than  the 
radiation  pulse  width,  and  the  effect  of  secondary 
avalanche  could  be  neglected  altogether. 

At  larger  E/p  when  the  ionization  -ate 
becones  larger,  the  discrepancies  in  the  ust  of  the 
complete  cross  .section  set  as  compared  to  the  limited 
set  are  still  significant.  Generally,  partiilos  that 
ire  not  programmed  to  undergo  momentum  translei 
interactions  rapidly  accelerate  to  velocities  where 
the  ionization  rate  Is  very  high.  However,  the 
particles  may  run  away  at  later  times  and  then 
contribute  little  to  Ionization,  Results  from  the 
complete  and  limited  cross  section  sets  show  this 
effect  in  Figure  8,  where  it  is  seen  that  the 
ionization  rates  are  large  enough  that  they  generally 


would  be  important  in  X-ray  pulse  experiments.  I'nis 
confirms  our  recommendation  that  momentum  transfer 
should  be  included  in  models  for  secondary  electrons, 
even  for  those  at  high  energies. 

The  results,  while  desired  for  should-be 

applicable  to  air  as  well.  The  de-exciiatlon  cross 
sections  for  Oj  partially  fill  in  the  depression  in 

the  energy  loss  rate  curve  in  Figure  2  between  the 
electron  energies  5  and  10  eV.  The  change  should  not 
influence  the  qualitative  nature  of  the  results 
presented  here.  Attachment  of  secondary  electrons  to 
the  electronegative  O2  molecule  is  an  effect  that  has 

been  neglected.  For  the  large  field  strengths  of 
Interest  in  AHSGEHP  it  would  not  appear  that  this 
effect  would  have  important  consequences  on  the 
results  presented  here.  There  is  no  intrinsic 
difficulty  in  incorporating  either  effect  into  the 
existing  algorithm. 


Figure  8.  Cumulative  Ionization  Rates  versus  Time  for 
the  Complete  and  limited  Cross  Section 

Sets.  Paraifieters  are  E  *  5  x  10^  V  x  m'^ 
and  p  «  1  torr.  500  histories  were 
followed  in  each  computer  run. 
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ATTN:  STRAT  &  THEATER  NUC  FOR  F  VAJDA 

DEPARTMENT  OF  THE  ARMY 

BDM  ADVANCED  TECHNOLOGY  CENTER 
AHN:  ATC-0  F  HOKE 
ATTN:  ATC-T 

FORT  HUACHUCA 

ATTN:  TECH  REF  DIV 


DEFENSE  ELECTRONIC  SUPPLY  CENTER 
ATTN:  DEFC-EAA 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  DB-6E 
ATTN:  DT-IB 
ATTN:  RTS-2B 

DEFENSE  LOGISTICS  AGENCY 

ATTN:  DLA-QEL  K  MASON 
ATTN:  DLA-SEE  F  HARRIS 

DEFENSE  NUCLEAR  AGENCY 

3  CYS  ATTN:  RAEV  TREE 

4  CYS  AHN:  STTI-CA 

DEFENSE  TECHNICAL  INFORMATION  CENTER 

12  CYS  ATTN:  DD 

FIELD  COMMAND  DNA  DET2 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  FC-1 

DNA  PACOM  LIAISON  OFFICE 
ATTN:  J  BARTLETT 

FIELD  COMMAND  DEFENSE  NUCLEAR  AGENCY 
ATTN:  FCPF  R  BLACKBURN 
ATTN:  FCPR 
ATTN;  FCTT 
ATTN:  FCTTWSUMMA 
ATTN:  FCTXE 

JOINT  CHIEFS  OF  STAFF 

ATTN:  CSS  EVAL  OFFICE  HDOO 

JOINT  DATA  SYSTEM  SUPPORT  CTR 
ATTN:  C-312R  MASON 
ATTN:  C-330 

JOINT  STRAT  TGT  PLANNING  STAFF 
ATTN:  JLKATTN:  DNA  REP 
ATTN:  JLKS 
ATTN:  JPPFD 
ATTN:  JPTM 


HARRY  DIAMOND  LABORATORIES 
ATTN:  CFAZi 
ATTN:  DELHD-NW-EAJMILEHA 
ATTN:  DELHD-NW-PTFLORY 
ATTN:  DELHD-NW-R  F  MCLEAN  22300 
2  CYS  ATTN:  DELHD-NW-RA'W  VAULT 
ATTN:  DELHD-NW-RH 
ATTN:  SCHLD-NW-P 
ATTN:  SLCHD-NW-EC 
ATTN:  SLCHD-NW-R 
ATTN:  SLCHD-NW-RA 
ATTN:  SLCHD-NW-RC 
ATTN:  T  TAYLOR 

NUCLEAR  EFFECTS  DIVISION 
ATTN:  R  WILLIAMS 

U  S  ARMY  ARMOR  &  ENGINEER  BOARD 
ATTN:  ATZK-AE-AR  J  DENNIS 

U  S  ARMY  BALLISTIC  RESEARCH  LAB 
ATTN:  SLCBR-VL  D  RIGOTTI 

USARMYBELVOIR  RD&ECTR 

ATTN:  DRDME-EJ  BOND  JR 

U  S  ARMY  COMMUNICATIONS  R&D  COMMAND 
ATTN:  DRSEL-NL-RO  R  BROWN 

U  S  ARMY  ELECTRONIC  TECH  DEV  LAB 
ATTN:  SLCET-SI  R  ZETO 

U  S  ARMY  ENGINEER  DIV  HUNTSVILLE 
ATTN:  HNDED-EDJ  HARPER 

U  S  ARMY  MATERIAL  TECHNOLOGY  LABORATORY 
ATTN:  DRXMR-BJ  HOFMANN 
ATTN:  DRXMR-HH  J  DIGNAM 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN:  LIBRARY 

US  ARMY  RESEARCH  OFFICE 
ATTN:  R  GRIFFITH 
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DEPARTMENT  OF  THE  ARMY  (CONTINUED) 

U  S  ARMY  STRATEGIC  DEFENSE  CMD 
ATTN:  BMDSC-AV  J  HARPER 
AHN:  BMDSC-HW  R  DEKALB 
AHN:  BMDSC-LEHRCWEBB 
ATTN:  DASD-H-SAV 

U  S  ARMY  STRATEGIC  DEFENSE  COMMAND 
AHN:  DASD-H-SA  R  C  WEBB 
AHM:  DASD-H-YD  F  HOKE 

U  S  ARMY  TEST  AND  EVALUATION  COMD 
ATTN:  AMSTE- 

U  S  ARMY  TRADOC  SYS  ANALYSIS  ACTVY 
ATTN:  ATAA-TFC  0  MILLER 

U  3  ARMY  WHITE  SANDS  MISSILE  RANGE 
ATTN:  STEWS-TE-AN  A  DE  LA  PAZ 
ATTN:  STEWS-TE-AN  J  MEASON 
ATTN:  STEWS-TE-AN  R  DUTCHOVER 
ATTN:  STEWS-TE-AN  R  HAYS 
ATTN:  STEWS-TE-N  K  CUMMINGS 
ATTN:  STEWS-TE-N  TARELLANES 
ATTN:  STEWS-TE-NT  M  SQUIRES 

USA  MISSILE  COMMAND 

ATTN:  AMSMI-LC-FS,  G.  THURLOW 
3  CYS  ATTN:  DOCUMENTS  SECTION  BLDG  4484 
ATTN:  DRCPM-PE-EAW  WAGNER 
ATTN:  HAWK  PROJ  OFCR  DRCPM-HAER 

XM-1  TANK  SYSTEM 

ATTN:  DRCPM-GCM-SW 

DEPARTMENT  OF  THE  NAVY 

LEAHY  CG  16 

ATTN:  WEAPONS  OFFICER 

NAVAL  AIR  SYSTEMS  COMMAND 
ATTN:  AIR  310 
ATTN:  AIR  350F 
ATTN:  AIR  54053 

NAVAL  AVIONICS  CENTER 

ATTN:  CODE  B415  D  REPASS 

NAVAL  INTELLIGENCE  SUPPORT  CTR 
ATTN:  NISC  LIBRARY 

NAVAL  OCEAN  SYSTEMS  CENTER 

ATTN:  CODE  9642B  TECH  LIB 

NAVAL  POSTGRADUATE  SCHOOL 

ATTN:  CODE  1424  LIBRARY 

NAVAL  RESEARCH  LABORATORY 

ATTN:  CODE  4040  J  BORIS 
ATTN:  CODE  4154  J  H  ADAMS 
ATTN:  CODE  4600  D  NAGEL 
ATTN:  CODE  4610  J  RITTER 


AHN:  CODE  4611  E  PETERSON 
ATTN:  CODE  4612  D  WALKER 
ATTN:  CODE  4612  RSTATLER 
AHN:  CODE  4613  A  B  CAMPBELL 
ATTN:  CODE  4614  L  AUGUST 
AHN:  CODE  4614  P  SHAPIRO 
ATTN:  CODE  4652  G  MUELLER 
AHN:  CODE  4653  A  NAMENSON 
ATTN:  CODE  4673  A  KNUDSON 
ATTN:  CODE  4682  C  DOZIER 
ATTN:  CODE  4682  D  BROWN 
ATTN:  CODE  5813  N  SAKS 
ATTN:  CODE  5313  W  JENKINS 
AHN:  CODE  5814  D  MCCARTHY 
AHN:  CODE  5814  M  PECKERAR 
ATTN:  CODE  5816  E  D  RICHMOND 
ATTN:  CODE  5816  R  HEVEY 
ATTN:  CODE  5816  R  LAMBERT 
AHN:  CODE  6810  J  KILLIANY 
AHN:  CODE  6816  H  HUGHES 

NAVAL  SEA  SYSTEMS  COMMAND 

ATTN:  CODE08K  NEWHOUSE 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  CODEF31 
ATTN:  CODE  F31  F  WARNOCK 
ATTN:  CODE  F31  K  CAUDLE 
AHN:  CODEH-20 
ATTN:  CODE  H23  R  SMITH 
ATTN:  F31 /DOWNS 

NAVAL  WEAPONS  CENTER 

ATTN:  CODE  343  FKA6A2,TECH  SVCS 

NAVAL  WEAPONS  EVALUATION  FACILITY 
ATTN:  CLASSIFIED  LIBRARY 

NAVAL  WEAPONS  SUPPORT  CENTER 
ATTN:  CODE  3073  T  ELLIS 
ATTN:  CODE  605  J  RAMSEY 
ATTN:  CODE  6054  D  PLATTETER 

NUCLEAR  WEAPONS  TNG  GROUP.  PACIFIC 
ATTN:  CODE  32 

OFC  OF  THE  DEP  ASST  SEC  OF  THE  NAVY 
ATTN:  LJABELLA 

OFC  OF  THE  DEPUTY  CHIEF  OF  NAVAL  OPS 
ATTN:  NOP985F 

OFFICE  OF  NAVAL  RESEARCH 

ATTN:  CODE  220  D  LEWIS 
ATTN:  CODE  414  L  COOPER 
ATTN:  CODE  427 

SPACE  &  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN:  CODE  5045.11  C  SUMAN 
ATTN:  CODE  50451 
ATTN:  NAVELEX  51024  C  WATKINS 
ATTN:  PME  117-21 
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DEPARTMENT  OF  THE  NAVY  (CONTINUED) 

STRATEGIC  SYSTEMS  PROGRAMS  (PM-1) 

AHN:  NSP-2301  M  MESEROLE 
ATTN:  NSP-2430J  STILLWELL 
ATTN:  NSP-2701 
AHN:  NSP-27331 
AHN:  NSP-27334 

DEPARTMENT  OF  THE  AIR  FORCE 

AERONAUTICAL  SYSTEMS  DIVISION,  AFSC 
AHN:  ASD/ENACC  R  FISH 
ATTN:  ASD/ENESS  P  MARTH 
AHN:  ASD/YH-EX  J  SUNKES 

AIR  FORCE  CTR  FOR  STUDIES  &  ANALYSIS 
2  CYS  ATTN:  AF/SAMI  TECH  INFO  DIV 

AiR  FORCE  GEOPHYSICS  l^BORATORY 
at;'N:  SULL 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 
ATTN:  AFIT/ENPC  BRIDGMAN 
ATTN:  LIBRARY 

AIR  FORCE  SYSTEMS  COMMAND 
ATTN:  DLCAM 
ATTN:  DLW 

AIR  FORCE  WEAPONS  UBORATORY,  AFSC 
ATTN:  NTAABCBAUM 
ATTN:  NTC  M  SCHNEIDER 
ATTN:  NTCASJ  FERRY 
ATTN:  NTCASJ  MULLIS 
ATTN:  NTCOXRTALLON 
ATTN:  NTCTGGOSS 
ATTN:  NTCTRKHUNT 
ATTN:  NTCTRRMAIER 
ATTN:  SUL 

AIR  FORCE  WRIGHT  AERONAUTICAL  LAB 
ATTN-  POE-2  J  WISE 
ATTN:  POOC-1  N  HARROLD 

AIR  FORCE  WRIGHT  AERONAUTICAL  LAB 
ATTN:  AFWAL/AADE 
ATTN:  AFWAL/MLPO  R  HICKMOTT 
ATTN:  AFWAL/MLTE 

AIR  LOGISTICS  COMMAND 

ATTN:  00-ALC/MMEDD/HARD  CONTROL 
ATTN:  OO-ALC/MMGR/OOC  CONTROL 

AIR  UNIVERSITY  LIBRARY 
atTn;  AUL-LSE 

BALLISTIC  MISSILE  OFFICE/DAA 
ATTN:  ENBE 
ATTN:  ENMG 
ATTN:  ENSN 

ATTN:  ENSN  M  WILLIAMS 
ATTN:  ENSN  W  WILSON 
ATTN:  SYST 
ATTN:  SYST  L  BRYANT 


BALLISTIC  MISSILE  0FFICE/A8RES 
AHN:  ENSN/HWARD 

ELECTRONIC  SYSTEMS  DIVISION 
ATTN:  INDC 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
ATTN:  TQTD  B  BALURD 

OFFICE  OF  SPACE  SYSTEMS 
AHN:  DIRECTOR 

OKLAHOMA  CITY  AIR  LOGISTICS  CTR 
ATTN:  DMM/R  WALLIS 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
AHN:  RBRJBRAUER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
ATTN:  ESRB  BUCHANAN 
ATTN:  ESRJ  SCHOTT 
ATTN:  ESRWSHEOD 

SPACE  DIVISION 

ATTN:  ALT 

SPACE  DIVISION 

ATTN:  YAS 

SPACE  DIVISION 

ATTN:  YAR  CAPT  STAPANIAN 

SPACE  DIVISION 
ATTN:  YD 

SPACE  DIVISION 
ATTN:  YE 

SPACE  DIVISION 
ATTN:  YG 

SPACE  DIVISION 
ATTN:  YN 

STRATEGIC  AIR  COMMAND 
ATTN:  INA 

STRATEGIC  AIR  COMMAND 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND 
ATTN:  XPFC 

STRATEGIC  AIR  COMMAND 
ATTN:  XPFS 

TACTICAL  AIR  COMMAND 
ATTN:  TAC/XPJ 

3416TH  TECHNICAL  TRAINING  SQUADRON  ATC 
ATTN;  TTV 

DEPARTMENT  OF  ENERGY 

DEPT  OF  ENERGY  ALBUQUERQUE  OPNS  OFC 
ATTN:  WSSB 
ATTN:  WSSB  R  SHAY 
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DEPARTMENT  OF  ENERGY  (CONTINUED) 

AEROJET  ELECTRO-SYSTEMS  CO 

AHN:  DTOOMB 

UNIVERSITY  OF  CALIFORNIA 

ATTN:  PLATHROP 

LAWRENCE  LIVERMORE  NATIONAL  LAB 

ATTN:  L-lOH  KRUGER 

AEROSPACE  CORP 

AHN:  L-13D  MEEKER 

ATTN:  D  FRESH 

1  ATTN:  L-156JYEE 

ATTN:  DSCHMUNK 

ATTN:  L-166  R  KALiBJiAN 

AHN:  G  GILLEY 

ATTN:  L-53  TECH  INFO  DEPT  LIB 

AHN:  IGARFUNKEL 

J  Af-rfc».  tirrkrsTJir* 

ni  1 IM.  TV  unvio 

ATT.N:  J  REINHEIMER 

i 

J  LOS  ALAMOS  NATIONAL  LABORATORY 

ATTN:  J  STOLL 

ATTN:  JWIESNER 

3  ATTN:  DLYNN 

ATTN:  J.B.  BLAKE 

i  ATTN:  MS  K  551  E  LEONARD 

AHN:  M  DAUGHERTY 

1  SANDIA  NATIONAL  UBORATORIES 

AHN:  PBUCHMAN 

AHN:  RSUUGHTER 

^  .  ATTN:  ORG  2100  B  L  GREGORY 

ATTN:  VJOSEPHSON 

S  ATTN:  ORG 21  low  DAWES 

ATTN:  W CRANE 

1  ATTN:  ORG  2115  J  EGOVER 

ATTN:  W  KOLASINSKI 

1  ATTN:  ORG  2143  TDELLON 

\  ATTN:  ORG  2144  P  V  DRESSENDORFER 

ALLIED  CORP 

j  ATTN:  ORG  2150  J  A  HOOD 

ATTN:  DOCUMENT  CONTROL 

1  AHN:  ORG  2320  J  H  RENKEN 

‘1  ATTN:  ORG  2321  L  POSEY 

ALLIED  CORP 

ATTN:  ORG  5143  J  L  DUNCAN 

ATTN:  EMEEDER 

i  AHN:  TWROBEL 

5  OTHER  GOVERNMENT 

AMPEX  CORP 

ATTN:  K  WRIGHT 

3 

[  CENTRAL  INTELLIGENCE  AGENCY 

\  ATTN:  OSWR/NEO 

AHN:  PPEYROT 

ANALYTIC  SERVICES,  INC  (ANSER) 

).  ATTN:  OSWR/STD/MTB 

ATTN:  ASHOSTAK 

J 

sj  DEPARTMENT  OF  TRANSPORTATION 

ATTN:  J  OSULLIVAN 

ATTN:  PSZYMANSK! 

J  ATTN:  ARD-350 

*i  NASA  GODDARD  SPACE  FLIGHT  CTR 

APPLIED  SYSTEMS  ENGRG  DIRECTOR 

ATTN:  JPRETZLERNUCSVMANG 

1  ATTN:  CODE  311.3  D  CLEVELAND 

1  ATTN:  CODE  313  V  DANCHENKO 

AVCO  SYSTEMS  DIVISION 

1  ATTN:  CODE  601  ESTASSINOPOULOS 

ATTN:  C  DAVIS 

1  ATTN:  CODE  660  JTRAI  NOR 

AHN:  D  SHRADER 

J  ATTN;  CODE  695  M  ACUNA 

ATTN:  F350DFANN 

1  ATTN:  CODE  701  W  N  REDISCH 

ATTN:  W  ERODING 

a  ATTN;  CODE  724.1  MJHABVALA 

i  NASA  LEWIS  RESEARCH  CTR 

BDM  CORP 

ATTN:  CMSTICKLEY 

1  ATTN:  MBADDOUR 

1  NASA  HEADQUARTERS 

BDM  CORP 

ATTN:  DWUNSCH 

1  ATTN:  CODE  DP  B  BERNSTEIN 

1  NATIONAL  BUREAU  OF  STANDARDS 

BEERS  ASSOCIATES,  INC 

ATTN:  B  BEERS 

1  ATTN:  C  WILSON 

2  CYS  ATTN:  H  W  BLOOMBERG 

1  ATTN:  CODE  A305  K  GALLOWAY 

1  ATTN:  CODE  A327  HSCHAFFr 

BOEING  CO 

0  ATTN-  CODE  A347J  MAYO-WELLS 

ATTN:  M/S  18-73  C  ROSENBERG 

ATTN;  CODE  A353  S  CHAPPELL 

ATTN:  M/S  2R-00  A  JOHNSTON 

ATTN;  CODE  A361 J  FRENCH 

ATTN:  M/S  2R-00  D  EGELKROUT 

1  ATTN:  CODE  C216J  HUMPHREYS 

AT-'N:  M/S  2R-00  EL  SMITH 

j  ATTN:  T  RUSSELL 

ATfN;  M/S2R-00  1  ARIMURA 

1 

j  DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ATTN:  M/S  2R-00  R  CALDWELL 

ATTN:  M/S  20-51  M  ANAYA 

i 

^  ADVANCED  RESEARCH  &  APPLICATIONS  CORP 

ATTN:  M/S  8J-78W  DOHERTY 

ATTN:  M/S  83-66  H  WICKLEIN 

1  ATTN:  L  PALKUTI 

ATTN:  M/S  88-43  OMULKEY 

1  ATTN:  RARMISTEAD 

ATTN:  M/S  9C-99  C  DIXON 

GENERAL  ELECTRIC  CO 

ATTN:  B  FLAHERTY 
ATTN:  G  BENDER 
ATTN:  LHAUGE 

GENERAL  ELECTRIC  CO 

ATTN:  CHEWISON 
ATTN:  nCOLE 

GENERAL  ELECTRIC  CO 
AHN:  D  PEPIN 

GENERAL  RESEARCH  CORP 
ATTN:  A  HUNT 

GEORGE  WASHINGTON  UNIVERSITY 
AHN:  A  FRIEDMAN 

GOODYEAR  AEROSPACE  CORP 

ATTN:  SEC  CONT  STATION 

GRUMMAN  AEROSPACE  CORP 
ATTN:  J  ROGERS 

GTE  GOVERNMENT  SYSTEMS  CORP 
AHN:  J  A  WALDRON 

HARRIS  CORP 

AHN:  EYOSt 
ATTN:  WABARE 

HONEYWELL,  INC 

AHN:  D  NIELSEN  MN  14-3015 
AHN:  RGUMM 

HONEYWELL,  INC 

AHN:  J  SCHAFER 
ATTN:  MS  725-5 

HUGHES  AIRCRAFT  CO 

ATTN:  B  CAMPBELL  El/EllO 

HUGHES  AIRCRAFT  CO 

ATTN;  RC  HENDERSON 

HUGHES  AIRCRAFT  CO 

ATTN:  WSCHENET 

HUGHES  AIRCRAFT  CO 

AHN:  MS-A2408J  HALL 

HUGHES  AIRCRAFT  INTL  SVC  CO 
ATTN:  ANAREVSKY 
ATTf5;  DSHUMAKE 
ATTN:  EKUBO 
ATTN:  WSC0TT 

IBM  CORP 

ATTN:  H  MATHERS 

IBM  CORP 

ATTN:  J  ZIEGLER 

IBM  CORP 

ATTN:  A  EDENFELD 
ATTN:  N HADDAD 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 

IIT  RESEARCH  INSTITUTE 
AHN:  IMINDEL 

ILLINOIS  COMPUTER  RESEARCH,  INC 
AHN:  ES  DAVIDSON 

INSTITUTE  FOR  DEFENSE  ANALYSES 
AHN:  TECH  INFO  SERVICES 

INTEL CORP 

AITN:  TMAY 

IRTCORP 

AHN:  J  AZAREWICZ 
AHN:  J  HARRITY 
AHN:  M  ROSE 
AHN:  MDC 
AHN:  N  RUOIE 
ATTN:  RMERTZ 

JAYCOR 

AHN:  MTREADAWAY 
ATTN:  R  STAHL 
AnN:  T  FUNAGAN 

JAYCOR 

AHN:  C  ROGERS 
AHN:  RPOLL 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  P  PARTRIDGE 
AHN:  R  MAURER 

JOHNS  HOPKINS  UNIVERSITY 

ATTN:  G  MASSON  DEPT  OF  ELEC  ENGR 

KAMAN  SCIENCES  CORP 
ATTN:  C  BAKER 

AHN:  DIR  SCIENCE  &  TECH  DIV 
ATTN:  J  ERSKINE 
ATTN:  N  BEAUCHAMP 
ATTN:  WRICH 

KAMAN  SCIENCES  CORP 
ATTN:  E CONRAD 

KAMAN  TEMPO 

ATTN:  DASIAC 
ATTN:  R RUTHERFORD 
ATTN:  W  MCNAMARA 

KAMAN  TEMPO 

ATTN;  DASIAC 

LITTON  SYSTEMS,  INC 

ATTN:  EL  ZIMMERMAN 
ATTN:  FMOTTER 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  rjUNGA  95-43 
ATTN:  J  SMITH 
ATTN:  REPORTS  LIBRARY 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  ABOROFSKY 
ATTN:  BKIMURA 


aHN:  EHESSEE 
AHN:  J  C  LEE 
AHN:  JCAYOT 
AnN:  L  ROSSI 
AHN:  PBENE 
AHN:  STAIMUTY 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
ATTN:  A.R;TOMME 
AHN:  LIBRARY 
AHN:  TECH  DATA  CENTER 

MAGNAVOX  ADVANCED  PRODUCTS  &  SYS  CO 
AHN:  WHAGEMEIER 

MARTIN  MARIEHA  CORP 
ATTN:  JTANKE 
AHN:  JIWARD 
ATTN:  MP-ieol  W  BRUCE 
AHN:  RGAYNOR 
ATTN:  TIC/MP-30 

MARTIN  MARIETTA  CORP 
ATTN:  T  DAVIS 

MARTIN  MARIEHA  DENVER  AEROSPACE 
AHN:  MS-b6074  M  POLZELLA 
AHN:  F'  KAk 
ATTN:  R  ANDERSON 
ATTN:  RESEARCH  LIBRARY 

UNIVERSITY  OF  MARYLAND 
AHN:  HCLIN 

MCDONNELL  DOUGLAS  CORP 
AHN:  APMUNIE 
AHN:  DLDOHM 
AHN:  M  STITCH 
ATTN:  R.L.  KLOSTER 

MCDONNELL  DOUGLAS  CORP 
AHN:  D  FITZGERALD 
ATTN:  P ALBRECHT 

MCDONNELL  DOUGLAS  CORP 
AHN:  TECH  LIB 

MESSENGER,  GEORGE  C 

ATTN:  G  MESSENGER 

MISSION  RESEARCH  CORP 
ATTN:  CLONGMIRE 

MISSION  RESEARCH  CORP 
ATTN:  R PEASE 

MISSION  RESEARCH  CORP 
ATTN:  J  LUBELL 
ATTN:  R CURRY 
ATTN:  WWARE 

MISSION  RESEARCH  CORP,  SAN  DIEGO 
ATTN:  J  RAYMOND 
ATTN:  V  VAN  LINT 
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DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 

MITRE  CORP 

AHN:  M  FITZGERALD 

MOTOROLA.  INC 

ATTN:  A  CHRISTENSEN 

MOTOROU,  INC 

AHN:  CLUND 
AHN:  LCURK 
AHN:  0  EDWARDS 

NATIONAL  SEMICONDUCTOR  CORP 
ATTN:  FC  JONES 

UNIVERSITY  OF  NEW  MEXICO 
AHN:  DRNEAMEN 

NORDEN  SYSTEMS.  INC 

AHN:  N  RIEDERMAN 
ATTN:  TECH  LIB 

NORTHROP  CORP 

AHN:  ABAHRAMAN 
AHN:  J  SROUR 
AHN:  ZSHANFIELD 

NORTHROP CORP 

ATTN:  EKING 
AHN:  SSTEV/ART 


RCA  CORP 

ATTN:  RSMELTZER 

RCA  CORP 

ATTN:  RKILLION 

RCA  CORP 

ATTN:  E  SCHMITT 
ATTN:  WALLEN 

RCA  CORP 

ATTN:  EVANKEUREN 

RENSSELAER  POL’/TECHNIC  INSTITUTE 
ATTN:  RGUTMANN 
ATTN:  RRYAN 

RESEARCH  TRIANGLE  INSTITUTE 
AHN:  M  SIMONS 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  AROVELL 
AHN:  GABO  TIC  LG  GREEN 
ATTN:  J  BELL 
ATTN:  J  BLIRSON 
ATTN:  KHULL 
ATTN:  VDE  MARTINO 
AiTN:  VSTRAHAN 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  TICAJOl 


PACIFIC-SIERRA  RESEARCH  CORP 

AHN:  H  BRODE  CHAIRMAN  SAGE 

PHYSICS  INTERNATIONAL  CO 
AHN:  J  SHEA 

POWER  CONVERSION  TECHNOLOGY  INC 
ATTN:  V  FARGO 

R  &  D  ASSOCIATES 

ATTN:  P  HAAS 
AHN:  W  KARZAS 

RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  P  DAVIS 
ATTN:  B BENNETT 

RAYTHEON  CO 

ATTN:  GJOSHI 
ATTN:  J  CICCIO 

RAYTHEON  CO 

ATTN:  A  VAN  DOREN 
ATTN:  HFLESCHER 

RCA  CORP 

ATTN:  G  BRUCKER 
ATTN:  VMANCINO 


ROCKWELL  INTERNATIONAL  CORP 
AHN:  LW  PINKSTON 
ATTN:  TIC  124-203 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  T YATES 
XhN:  TICBA08 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DLONG 
ATTN:  DMILLWARD 
ATTN:  DSTROBEL 
ATTN:  J  NA8ER 
ATTN:  JSPRATT 
AHN:  L  SCOTT 
ATTN:  RFITZWILSON 
AHN:  R  J  BEYSTER 
ATTN:  V  ORPHAN 
ATTN:  VVERBINSKI 

SCIENCE  APPLICATIONS  INTL  COPP 
ATTN:  .W  CHAOSEY 

SCIENTIFIC  RESEARCH  ASSOC.  INC 
ATTN:  HGRUBIN 

SPERRY  CORP 

ATTN:  J INDA 
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DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 


SPERRY  CORP 

AHN:  PMARROFFINO 

SPERRY  CORP 

AHN:  DSCHOW 

SRI  INTERNATIONAL 

ATTN:  APADGEH 

SUNDSTRAND  CORP 

AHN:  C  WHITE 

SYSTEM  DEVELOPMENT  CORP 

AnN:  PRODUCT  EVAL  LA3 

SYSTRON-DONNER  CORP 
AHN:  J  RAY 

TELEDYNE  SYSTEMS  CO 
ATTN:  RSUHRKE 

TEXAS  INSTRUMENTS,  INC 

ATTN:  E  JEFFREY  MS  961 
AHN:  FPOBLENZMS3143 
ATTN:  J  SAL2MAN 
AHN:  R  MCGRATH 
AHN:  T  CHEEK  MS  3143 


TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  AWinELESMSRl/2144 
ATTN:  D  CLEMENT 
AHN:  FFRIEDT 
ATTN:  H  HOLLOWAY 
AHN:  MS  ASH 
2CYS  AHN:  0  ADAMS 
ATTN:  PGUILFOYLE 
AHN:  PR  REID 

2CYS  AHN;  R  PLEBUCH  HARDNESS  &SURV  UB 
ATTN:  RVOMHAHEN 
AHN:  TECH  INFO  CTR 
AHN:  W  B  ADELMAN 
ATTN:  W  ROWAN 
AHN:  W  WILLIS 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
AHN:  CSLASNEK 
AHN:  FFAY 
ATTN:  J  (50RMAN 

WESTINGHOUSE  ELECTRIC  CORP 
AFTN:  D  GRIMES 
AHN:  HKALAPACA 
AHN;  RCRICCHI 

WESTINGHOUSE  ELECTRIC  CORP 
ATTN:  SWOOD 
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